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Abstract. We derive a method to calculate the multi-channel K matrix applicable to a broad class of 
models in which mesons linearly couple to the quark core. The method is used to calculate pion scattering 
amplitudes in the energy region of low-lying Pll and P33 resonances. A good agreement with experimental 
data is achieved if in addition to the elastic channel we include the irA and aN (o"Z\) channels where the 
cr-meson models the correlated two-pion decay. We solve the integral equation for the K matrix in the 
approximation of separable kernels; it yields a sizable increase of the widths of the zi(1232) and the 
Af(1440) resonances compared to the bare quark values. 
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1 Introduction 

Among the low-lying nucleon excitations, the Pll Roper 
resonance 7V(1440) as well as its counterpart in the P33 
partial wave, the Z\(1600), play a special role due to their 
relatively low masses and due to the rather peculiar be- 
havior of the scattering and electro-excitation amplitudes. 
The constituent quark model (CQM) in which the excited 
' states are treated as bound states yields partial decay 
widths which are generally considerably smaller than the 
experimental values [1-8] unless all widths are scaled by a 
factor of 2-3 in order to fit the experimental width of the 
Zi(1232). This indicates that the structure of the Roper 
■ can not be explained by a simple excitation of the quark 
' core (like most of the other low-lying states) and that 
other degrees of freedom need to be included. 

Several attempts to understand and explain the na- 
ture of the Roper resonance have been proposed and dis- 
cussed. They include the study of the Roper resonance on 
the lattice [9-12], investigations in different quark mod- 
els, like the models based on qqqqq configurations [13- 
15] or those incorporating the meson cloud consisting pri- 
marily of pions and cr-mesons [16-23], as well as in mod- 
els with hybrid (qqqg or glueball) configurations [24,25]. 
Pion-nucleon scattering in the region of the Roper reso- 
nance has also been studied in the framework of chiral 
perturbation theory [26-28]. The scattering amplitudes 
and the decay rates have been well established in cer- 
tain phenomenological approaches, for example in energy- 
independent partial-wave analysis of ttN — > t:t:N scat- 
tering [29-31], energy-dependent analysis [32,33], model- 
independent analysis [34] , as well as in dynamical coupled- 
channel models [35,36] and in effective-Lagrangian mod- 
els (see, for example, [37-44]). Such models typically con- 



tain a large number of ingredients (e.g. coupling constants, 
cut-offs, meson and baryon states) in order to parameter- 
ize the observed scattering (or electro-production) ampli- 
tudes and to simultaneously fit many measured observ- 
ables. However, because of the multitude of model ingre- 
dients, it is usually difficult to determine which degrees of 
freedom are truly relevant for a particular resonance and 
to obtain insight into its quark structure. Furthermore, 
a similar level of agreement with data can be achieved 
using rather disparate sets of parameters; using results of 
quark-model calculations may offer an important guidance 
in particular to the choice of the form factors of resonant 
states. 

The first aim of the present work is to develop a method 
which would allow solutions obtained in a broad class 
of quark-model calculations using bound-state boundary 
conditions to be incorporated into a dynamical framework 
that could be used to predict observables measured in me- 
son scattering experiments as well as in electro- weak pro- 
duction of mesons on the nucleon. The interplay of dif- 
ferent baryonic and mesonic degrees of freedom in such 
processes can in many cases considerably alter the results 
obtained by simply calculating the transition form-factors 
of the baryons in the underlying quark model. 

In addition, our goal is to set up a computational 
scheme beyond the simple approximation for the K matrix 
and connect the observables extracted in our analysis of 
the relevant degrees of freedom to the predictions of the 
underlying quark model, which would eventually give us 
a chance to discriminate between different models. 

Our method is not intended to be a competitive ap- 
proach to those using effective Lagrangians. Rather, the 
method complements them by establishing a link between 
the effective models and the underlying baryon structure. 
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Its main distinctions with respect to the (non-dynamical) 
effective Lagrangian approaches are: (i) baryons are treated 
as composite particles from the very beginning; the strong 
and electro-weak form-factors are derived from baryon in- 
ternal structure and not inserted a posteriori; as a conse- 
quence the method introduces a much smaller number of 
free parameters; (ii) the meson cloud around baryons is 
included in a consistent way also in the asymptotic states; 
(iii) through the Kohn variational principle for the K ma- 
trix it allows one to determine the internal structure of the 
baryon resonances which may substantially differ from the 
structure determined in the calculations using bound-state 
boundary conditions. On the other hand, because of the 
complex structure of the baryons in the model, we are not 
able to treat a very large set of ingredients. Consequently, 
the approach in its present form cannot be used at higher 
energies where many channels open. 

In several aspects the present method is similar to 
the method of solving the integral equation for the T 
matrix described in [35] and used in [36]. In both ap- 
proaches the coupled-channel equations are solved beyond 
the usual Born approximation including the effect of the 
meson cloud in the baryon; as a consequence, in both cases 
the expression of the T matrix (K matrix) involves the 
physical nucleon rather than a structureless particle of 
effective Lagrangian approaches. The starting points are 
however different: the approach of [35] starts from baryon- 
meson and meson-meson Lagrangians and using the uni- 
tary transformation method derives a model Hamiltonian 
with dynamically generated baryon-meson vertices which 
can be related to the vertices calculated in quark models. 
Our approach is complementary in that it starts from the 
quark-model type calculation and extends the many-body 
wave-function such that it properly includes the asymp- 
totic boundary conditions of a particular channel. The def- 
ficicncy of our approach in the present form is the lack of 
meson self-interactions. 

Our future plan is to extend the approach to include 
meson photo- and electro-production which requires the 
inclusion of new channels. In order not to deal with an 
excessively large number of ingredients in that case, the 
aim of the present analysis is also to reveal which degrees 
of freedom are most relevant for an explanation of the 
structure and the dynamics of the low-lying resonances 
in the Pll and P33 partial waves and to pin down the 
parameters governing these resonances. 

In our previous work [45] we have developed a method 
to calculate the scattering and the pion electro-production 
amplitudes in models in which the pions linearly couple 
to the quark core. In such models it is possible to find the 
exact expression for the T matrix as well as for the K ma- 
trix without explicitly specifying the form of the asymp- 
totic states. The method has been successfully applied to 
the calculation of the phase shift and the pion electro- 
production amplitudes for the P33 partial wave by consid- 
ering only the elastic channel. In this paper we extend the 
method to calculate directly the multi-channel K-matrix. 
The resulting matrix is symmetric and real, thus ensur- 
ing unitarity. Since we are now considering processes at 



relatively higher energies we generalize the approach by 
adopting fully relativistic kinematics. 

The first attempt to treat dynamically the resonant 
pion scattering in a quark model with excitations of the 
quark core has been done in the framework of the Cloudy 
Bag Model (CBM). For the P33 partial wave [19,23] the 
experimental phase shift in the energy range of the 4(1232) 
resonance has been well reproduced. In the 7V(1440) case 
the width of the resonance has been considerably under- 
estimated [17,19,22] except in the case of the coupled- 
channel calculation [21] using the K-matrix approach in- 
volving the 7rZ\ and tt- Roper inelastic channels in the static 
approximation. The latter approach is similar to the me- 
thod of our sect. 3.3 except that we use relativistic kine- 
matics. 

In the present work we restrict ourselves to a class of 
Hamiltonians with linear meson-baryon coupling, neglect- 
ing meson self-interaction. We are particularly interested 
in the interplay of the degrees of freedom involving the 
A isobar and the cr-meson which we expect to play the 
dominant role in the dynamics of the Roper resonances; 
in this respect our work is similar to that of Rcf. [39, 
46]. We were, however, not able to confirm the most in- 
triguing conclusion found in [39] that the formation of the 
iV(1440) resonance can be explained without introducing 
a bare Roper state. 

In sect. 2 we introduce a general form of models in 
which mesons linearly couple to the quark core. The in- 
put to our computational scheme are the matrix elements 
of the meson interaction between the quark states. We 
construct the K matrix and the corresponding principal- 
value states which incorporate in a systematic way the 
many-quark quasi-bound states describing the excitations 
of the quark core. We assume that the decay into two pi- 
ons which dominates the inelastic processes in the lower 
energy region proceeds through a two-body decay involv- 
ing either an unstable isobar which in turn decays into the 
nucleon and the second pion or an unstable meson and the 
nucleon. 

In sect. 3 wc show how a quasi-bound state calculated 
in the underlying quark model can be inserted in a state 
that obeys proper scattering boundary conditions. We de- 
rive a set of coupled integral equations for the meson am- 
plitudes and the parameters governing the strength of the 
quasi-bound states which are responsible for the resonant 
parts of the scattering amplitudes. We first solve the set in 
a simple approximation equivalent to the so-called Born 
approximation for the K matrix widely used in various 
multi-channel approaches. In order to find the solution 
beyond the Born approximation we introduce an approx- 
imation which makes the kernels separable, however, en- 
suring that they reduce to the exact form when evaluated 
on-shell, as well as preserving the proper symmetries of 
the K matrix. 

The results of our analysis using the quark wave-func- 
tions from the Cloudy Bag Model are presented in sect. 4. 
We analyze to what extent various approximations are 
able to explain the resonant behavior of the Pll and P33 
partial amplitudes. We stress the important role of the 
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background terms originating from the crossed (u-channcl) 
processes, and show that the inclusion of the aN channel 
considerably improves the results in the Pll case. The 
meson-baryon coupling constants needed to fit the scat- 
tering amplitudes in the Born approximation for the K 
matrix are generally considerably larger than those pre- 
dicted by the quark model. Solving the integral equation 
for the K matrix beyond the Born approximation we show 
that the resulting dressing of the vertices as well as the 
influence of the neighboring resonances can explain the 
strong enhancement of the coupling constants. 



2 Coupled-channel K-matrix formalism 

2.1 The K matrix for models with linear meson-quark 
coupling 

We consider here a rather broad class of models in which 
the mesons linearly couple to the three-quark core. The 
part of the Hamiltonian referring to mesons can be written 

as 

H n = dk^2{^kaj mt (k)ai mt (k) 

J Imt 

+ [Vi mt (k)a lmt (k) + Vl t {k)a\ mt {k)] } , (1) 



where a\ mt (k) is the creation operator for a ^-wave meson 
with the third component of spin m and - in the case of 
isovector mesons - the third component of isospin t. In 
the case of the p-wave pions, the source can be cast in the 
form 

3 



V mt (k) = -v(k)J2<rt 



(2) 



?=i 



with v(k) depending on a particular quark model and con- 
taining the information about the underlying quark struc- 
ture. (We omit the index I from the pion operators.) We 
assume that the interaction V(k) can generate bare quark 
states with quantum numbers different from the ground 
state by flipping the spin and isospin of the quarks, and 
furthermore, excite quarks to higher spatial states. In par- 
ticular, the state with the flipped spins and isospin cor- 
responds to the 4(1232) isobar, while exciting one quark 
from the Is to the 2s state generates an excited three- 
quark state associated with the Roper resonances iV(1440) 
or 4(1600). 

Chew and Low [48] considered a model similar to (1) 
except that they did not allow for excitations of the nu- 
cleon core. They showed that the T matrix for 7riV scat- 
tering was proportional to ('P^ (W)\V m t(k)\ l I / N} , where 
&(~\W) were the incoming states. In general, the corre- 
sponding formula for the K matrix cannot be written in 
such a simple form. However, in the JT basis, in which the 
K and T matrices are diagonal, it is possible to express 
the K matrix as 



kW 



W T {W)\\V{kW N ). (3) 



The corresponding principal- value (PV) state obeys a sim- 
ilar equation as the in- and out-going states in the Chew- 
Low model and takes the form: 



W?t(W))=M { [a\k )\V N )] 
V 



JT 



[V(k )\V N )] 



JT 



H-W 

where [ } JT denotes coupling to good J and T, and 



uj E_ 



N 



k W 



(4) 



(5) 



We work in the center-of-mass system, W is the invariant 
energy of the system, ujq and fco 



m 2 are the 

energy and momentum of the pion, and En is the nucleon 
energy: 



wo 



W -E N 



W 2 



M N + ml 



2W 



E 



N 



^M N + kl. 



It is worthwhile to notice that a PV state is a su- 
perposition of incoming and outgoing waves {i.e. can be 
regarded as a standing wave); the in- and out-states be- 
come meaningful only when the T and the S matrix are 
constructed from the K matrix: 



T = 



K 



(1-iK)' 



S=l + 2iT. 



(6) 



2.2 Channels including the pion and unstable isobars 

We now extend the above formulas to the multichannel 
case. We apply the usual approach used in phenomeno- 
logical analyses (e.g. [29]) in which it is assumed that the 
two-pion decay proceeds through some intermediate un- 
stable particle, either a meson or a baryon. 

Let us first consider the situation in which the inter- 
mediate particle is a baryon, B, (e.g. the 4 isobar which 
dominates at low energies) which in turn decays into the 
nucleon and the second pion. The decays into the nucleon 
and an unstable meson will be treated in the next subsec- 
tion. 

In analogy with (4) we introduce the PV state corre- 
sponding to the ttB channel as 



\$fr{W,M)) =M<j [a\kx)\$ B {M)) 
V 



JT 



H-W 



V(fci)|«MM)> 



JT 



(7) 



Here k\ is the momentum of the first pion; the energy E 
of the isobar B and the energy of the first pion are related 
through 



E 



W 2 + M 2 - rnj 
ZW 



uj 1 =W-E, 



(8) 
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where M is the mass of B, i.e. the invariant mass of the 
nucleon and the second pion, and 



M 



u>\E 
~k~iW ' 



(9) 



The normalization of the PV states (4) and (7) states 
is dictated by their respective first terms which represent 
a free pion and the nucleon or the B isobar. The states 
\P a , a — N,(B, M), . . . are then normalized as 

(V a (W)\¥\W>)) = 8(W - W>) [5 a ^ + K 2 a ,p] 

(see e.g. [47], eq. (7.28)). In the case of the 7riV channel 
a = N is a single index, while in the wB channel it includes 
also the invariant mass M, so S a . a should be interpreted 
as S(M — M'). The PV states are not orthonormal; the 
orthonormal states are constructed by inverting the norm: 



\W a {W)) =J2 [l + K 2 



-1/2 



fia 



\* P (W)) (10) 



(see e.g. [47], eq. (7.29)). The state \P B (M) in (7) repre- 
senting the intermediate B baryon is constructed using 
(10), and is therefore normalized as (§ B (M)\& B (M')) = 
6(M — M'). The construction of the orthonormal states 
corresponding to intermediate isobars is discussed in Ap- 
pendix B. 

For a process in which the initial pion-nucleon sys- 
tem with invariant mass W decays into the (first) pion 
with momentum k and the pion-nucleon system with the 
quantum numbers of the intermediate baryon B, we write 



mesons is straightforward but more complicated because 
of a larger number of involved channels. 

The (T-meson appears as a chiral partner of the pion in 
several versions of chiral quark models which also provide 
the form of its coupling to the quark core, e.g. in the lin- 
ear cr-model with quarks [50-52] or in various bosonized 
versions of the Nambu-Jona-Lasino model [53-55] . In the 
non-linear realizations of the models its coupling is real- 
ized by two correlated pions, see e.g. [18,37]. 

In this work we use a purely phenomenological ap- 
proach and take the form (1) for the cr-meson coupling. 
We assume that only s-wave cr-mesons couple to the quark 
core such that the interaction vertex takes the form: 



V lt (k) = V»(k)w <r (n), V»(k)=G a 



k 



Apart from the momentum k, the one-cr-meson states are 
labeled by the invariant mass of the two-pion system, 
2m n < yu < oo. Here uj 2 k = k 2 + [i 2 and w a (\i) is a nor- 
malized mass distribution function centered around the 
nominal value of the cr-meson mass with the correspond- 
ing width, modeling the resonant decay into two pions. 

In analogy with (3), (11), (12), and (13), we first intro- 
duce the matrix elements of the K matrix referring to the 
aN channel in the case of the Pll partial waves. They in- 
volve the cr-meson with the energy and momentum to^ and 

fc M and the nucleon with the energy E 1 ^ = + k^ on 

the one side, and on the other side either the pion (loo, k ) 
and the nucleon, the pion (uj\, ki) and the intermediate 
baryon with the invariant mass M, or another cr-meson 
and the nucleon: 



K% wN (k,k ,M) = -n^^(*?AW)\\V(kWB(M)). 

(11) 

The irB to irN transition matrix element of the multi- 
channel K matrix is 



K^ B (k, ki, M) = -^^{wfriW, M)\\V{k)\\W N ) . 

(12) 

The matrix element corresponding to the tiB to ttB' tran- 
sition (e.g. Tr(fci) + A(M) -> n(k) + A(M' j) is given by 

K^, nB (k,h,M',M) 

= -*sJ^^fT(W,M)\\V(kW B ,(M')). (13) 



2.3 Channels including the cr-meson 

We consider here the situation in which the decay proceeds 
through an unstable meson and a baryon. We derive the 
expression for the K matrix in the simplest case of the 
aN channel (in the Pll wave) and the Art channel (in 
the P33 wave) which dominate the inelastic processes in 
the energy range below ~ 1700 MeV for these two par- 
tial waves. The extension of the method to other unstable 



K^ N (W^) = -nM, <!2fi(W)|^(MI<M , 



K^ aN (W,u) = -ttM, mAW,»)\\V(kv)\\V N ) , 



2 2 

kILb(W,»,M) = -tt^ (^? h (W,M)\W(k^ N ) , 



K^ aN (W,M,n) = -ttM (n i (W, f j)\\V(k 1 )\\# B (M)) , 



K 



aNaN 



(W,n,n') = -^(nx(W,ii')\V^(k^ N ) , 



(14) 



where TVq and N\ have been defined in (5) and (9), 



and 



UJ,, 



W 2 



M 2 N + ^ 2 



2W 



w 



(15) 



In the P33 partial wave the decay into two correlated 
I = pions can proceed only through the intermediate s- 
wave cr-meson and the A isobar. In this case the channel is 
labeled by the invariant mass of the two correlated pions, 
/i, and the invariant mass of the pion-nucleon system from 
the A resonance, M. The corresponding elements of the 
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K matrix contain four invariant masses which is rather 
difficult to handle computationally. In sect. 3.4 we discuss 
a method which simplifies the calculation by averaging 
over the A invariant mass such that the dependence on 
M is eliminated. The matrix elements then assume the 
same form as (14) with N replaced by A, Mm by M, and 



E N by W = ^M 2 + fe2. 



In the PI I and P33 case the p-meson can be included in 
a similar way but since the experimental data [56] indicate 
that the contribution of the p-meson is almost negligible 
in the energy region under consideration we do not include 
it in our calculation. 



2.4 Constructing the T and the S matrix 

The multichannel K matrix acquires the following form 



Knn 


Knb(M') 




K BN (M) 


K BB ,(M, M') 


K Ba (M^>) 




K aB {n,M) 


K aa {ji,ii') 



where we have used a shorthand notation B for the ttB 
channel and a for either the aN or the a A channel. The 
T matrix is related to the K matrix through (6). In the in- 
elastic channel the matrix elements depend on the contin- 
uous variable (M or fi), yielding a set of coupled integral 
(Heitler) equations: 

T NN = K NN + iT NN K NN 

nW—rrin 

+ dMT NB (M)K BN (M) 

rW-M N 

+ i / dfj,T N „(fj,)K aN (n) , 

J2m 7r 

T NB {M) = K NB {M) + \T NN K NB {M) 

_ rW—m-K 

+ &M'T NB ,{M')K B , B {M',M) 

r-W~M N 

> I d//7 A >(//)A'„/ ; ;/'..\/;. 



2m, 

T N <j(li) = K Na (n) + iT NN K Na (fi) 

+ i£ / dMT NB (M)K Ba (M,ft) 

pW—M N 

+ i d^ 1 ' T N<r (n')K air (n' » . 

J2m 7r 



(16) 



In the P33 case the nucleon mass in the integral over 
fj, is replaced by the W-^-dependent averaged invariant mass 
of the intermediate A. 

The unitarity of the S matrix it fulfilled provided the 
K matrix is real and symmetric which is especially impor- 
tant when we use approximate methods; in such a case it 
is considerably more advantageous to use a certain ansatz 
for the K matrix (or, equivalently for the principal-value 
state) rather than for the T matrix since in the latter case 
the unitarity has to be enforced at each step of the calcula- 
tion. Let us remark that for a general chiral quark model, 



the K matrix and the corresponding principal-value state 
can be calculated variationally using the Kohn variational 
principle 

(5V P \H- W\V P ) =0, (17) 
where & v is a suitably chosen trial state. 

3 The integral equations for the K matrix 

3.1 Ansaetze for the channel PV states 

In the formal expressions for the PV states (4) and (7) the 
interaction V(k) generates bare quark states with quan- 
tum numbers different from the ground state, as well as su- 
perpositions of bare quark states and one or more mesons. 
We choose a particular ansatz which implies the proper 
asymptotic behavior of different channels consisting of a 
meson and a baryon carrying its own meson cloud. The 
form of the pion state in a ttB' channel can be read-off 
from the general relation (55) holding for the eigenstate 
of the Hamiltonian (1). Multiplying (55) for \P A = &f T by 

(\P B i(k)\ we obtain 



(w fc +E B ,{k) 



W)(V B ,{k)\a mt {k)\Vf T {W)) 
-{* B >{h)\vUWfrm) , (18) 



with E B i(k) — \/M' 2 + k 2 , where M' is either the mass 
of a stable baryon or the invariant mass of the ttN sub- 
system corresponding to the intermediate baryon B' . For 
the elastic channel we assume the following ansatz: 

\*jt(W)) - Afo {[a t (fco)|<Mfco))]' 7T + £ <%(W)\<hi) 



+ / dk 



I 

+J2 J dM J dk 



« t (fc)l'Mfc)>] 



Lu k + E N (k)-W 



JT 



+ / dfi / dk 



uj k + E B (k)-W 
X3r(fc,*b,/i) 



[at(fc)|!MM)>] 



JT 



Lu^+Ejrik) - W 



(19) 



The first term, as discussed in the previous section, defines 
the channel and determines the normalization, the second 
term is the sum over bare quark states, denoted as 
with quantum numbers of the channel. The amplitudes 
X are proportional to the amplitudes (18); the first term 
corresponds to the one-pion state on top of the ground 
state, the terms in the sum to the one-pion states around 
different excited states, and the last term to the one-cr 
state around cither the nucleon (Pll) or the A (P33), 
with Ejt (k) denoting the energy of the respective baryon. 
Above the one- and two-pion thresholds these amplitudes 
are related to the elastic and inelastic elements of the on- 
shell K matrix: 



K NN {W) ^t:M^ X jt (to,**)), 
K BN (W,M) = TtN N lX Bh 
K aN { W, n) = 7T A/uA^ X jt , ko, n) ■ 



JT (fci,fc ,M), 



(20) 
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In the Pll case, one of the <Pn states is the nucleon. 

The ansatz can be simplified by realizing that the main 
contribution to the integrals over the invariant mass M, in 
(19), comes from M close to Mb, i.e. the resonant energy 
of the isobar B. In appendix B we show that in such a 
case one can write the state &b m a simplified form 



\*b{M))=wb{M)\# b {M)), 



(21) 



where SP B is dominated by the bare quark configuration 
and only weakly depends on M. The function wb(M) can 
be identified with the mass distribution function a(M) 
multiplied by the kinematic corrections - the Blatt-Weiss- 
kopf barrier-penetration factor - that ensure proper thresh- 
old behavior (see e.g. [29]). Using this approximation the 
integration over M in (19) selects M = Mb and similarly 
fi = m a in the last term, yielding 



n 



\*jt{W)) = M \[a\k )\* N {k ))] JT + Y j <%{W)\$k) 

x3t (*>*o) 



+ / dfc 



/ 

B 

I 



u k + E N {k)-W 
X^(k,k ,M B ) 



dk 



+ / dfc- 



u k + E B (k)-W 
XjT(k,k ,m, 



aHk)\4r N (k))] JT 

a\k)\$ B {M B ))] JT 
b\k)\$ JT )\, (22) 



uj k +E JT {k)-W 
where ui k — \Jk 2 + m 2 a . We have introduced 

x ™ = w B (M)x^ and *Sr = «v(m)x3t ■ 

The inelastic channels corresponding to the irB and 
aN channels can be written in compact forms 

\*f T (W,M)) 

= Af B w B {M){[a){k 1 )\$ B {M))] JT + ^c^(W, M)|*r) 
, j dk X^(k,h,M) 



W(k)\# N (k))] 



uj k + E N (k) - W 

B' 



Lu k + E B >(k) - W 



+ 



(23) 



and 

\*5t{W,ii)) 

Xjt (^i A 4 ) 



n 



+ / dfc 



w/c + S w (fc) 



[at(fc)|^(fc))] 



/ 

X^(fc,fc^,m CT ,/i) 



dfc 



£ fe + £j T (fc)- W 



(24) 



Here we have assumed the following factorization 



Xjt^ = w H {m H )w H ,{m' H )x 



1 \cyH'H 
JT > 



(25) 

where iJ stands for cither the ttB channels or aN, and 
niH is the invariant mass (either M or y(i) . Above the two- 
pion threshold the meson amplitudes are related to the K 
matrix by 

K H , H (W,m' H ,m H ) = itNh'Nh Xjt" (kn' , k H ,m' H ,m H ) . 

The requirement that the K matrix be symmetric imposes 
the constraint Xjt H (k' > k) = Xjt (k, k') on the pion am- 
plitudes. The ansaetze therefore ensure that the scattering 
amplitudes are directly proportional to the on-shell K ma- 
trix needed in the equation for the T matrix (16). For the 
off-shell matrices, equality (26) is not fulfilled in general. 



3.2 Derivation of the coupled integral equations 

The equations for the pion amplitudes x an d the coeffi- 
cients c-jz in the ansaetze (19), (23) and (24) are derived ei- 
ther from (17) or directly from the commutation relations 
(55) which hold for our particular choice of the quark-pion 
interaction. 

By requiring stationarity with respect to the variation 
of the coefficients c^>, and we get 



(W 



M° n )c%(W) 
= V Nn (ko) - 



dk 



u k + E N (k)-W 



+ 



+ 



i 



dk 



V^(k)x^T N (k,ko,M BI ) 



dfc^ 



u k + E B >(k)-W 
{k)x a jT(k,kQ,m a ) 



u k + E JT (k)-W 



(27) 



(W-M&c&WM) 



M 
BIZ 



(ki) + J 



dk 



+ 



R' J 



dk 

+ J dk VjreVI>MjT 



^vR(fc)xy T B (fc,fci,M) 

w k + E N (k)-W 



uj k + E B '{k)-W 



ti k + E JT (k)-W 



(28) 



{W-M^)^{W,n) 



+ 



R' J 



J 



+ / dfc 



uj k + E N (k)- W 

dh V^n(k)x^(k,k M M B ',fi) 
uj k + E B >(k)-W 

Kn(k) x7 T {k,k^,m a ,^) 
Co k + E JT {k)-W 



(29) 
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Here: 



V NU {k) = {$n\\V{k)\\V N ) = Z- 1/2 {$ n \\V(k)\\$ N ) , 
V^(fc) = (4>n\\V(k)\\$ B (M)} = Z B 1/2 (<P K \\V(k)\\<P B ) , 

v£ w (fc) = {<t>n\\v»{k)\\$ N ) - ^ 1/2 <<%||v^(fc)||<M , 

where Z B is the wave-function normalization, while 
($n\\V(k)\\$ N ) and (<P tc ||V(A;)||#b> are obtained from 
the underlying quark model. In the P33 case Vfi n (k) is 
replaced by V% n (k). 

In the Pll case, one of the states <Pti is replaced by 
the (exact) ground state in which the requirement of sta- 
tionarity is equivalent to the requirement that the channel 
states are orthogonal to the ground state for W > Mm- 
In this case the mass of the bare state in (27), (28) 
and (29) is replaced by the ground-state mass Mjv, while 
the matrix elements are given by 



canceling wb(M)wb(M') on both sides, 

X% B (k,k l7 M\M) 

= K B M P M {kM) J2 ? n(W,M)Vif n (k) 



V Nn (k) -> (W-M N ) 



Vik(k) - (W-M N ) 



($ N \\V(k)\\V N ) 
Lu k + E N (k)-W 

{* N \\V(k)\\$ B (M)) 
co k + E B (k)-W 



(30) 
(31) 



Requiring stationarity with respect to pion amplitudes 
leads to the familiar Lippmann-Schwinger equation for the 
K matrix. The equation for the Xjt amplitude which is 
related to the elastic part of the K matrix is obtained 
from (18) for B = B' = N . Using our ansatz (22) and 
taking Mj = Mt = \, we obtain, after multiplying (18) 

by Ci !_ TOlTO Ci I _ tlt and summing over m and t, 



n 

,lC NN (k,k')x^(k>,k ) 



dk 



Lo' k + E N (k')-W 
,IC^ B B (k,k')x^(k',k ,MB) 



u> k + E B (k<)-W 
,!C^l{k,k')x^{k'M,m a ) 



where we have introduced the kernels 
K N M B {k,k') = - (#N(k)\al, tl (k') 

mtm't' 

vUk) + K + E N (k) - W)a mt (k)] \$ B (M)) 

j 1 7 1 1 71 7^1 

XU J B i-m'lm' L 'T B i-t'U' L 'ii-mlm ii-tlt • V JJ J 

For B = N, $ B {M) reduces to <P N and M to M N . 

For the general amplitude involving the irB channels 
we use our ansatz (23) which yields, after using (25) and 



B" J 

+ I ik - 



n 

B'B" 

M'Mrj 



„(fc,fc / )xfr B (fc , ,Ai,M B »,M) 



u' k + E B »(k')-W 



Here 



K 



mm' (k, k') 



Cj' k +E JT {k>)~W 



]T ($B{M)\al, t ,{k') 



(34) 



Vl t (k) + (uj k + E{k) - W)a mt {k) \9 B , (M')) 



v 2 /~< J 2 (~< J 2 r<^ 2 

A u J^-m'lm' i -fit' ° J B i -mlm u T B I -tit ■ 

The form of (34) justifies the factorization (25) for the 
amplitude Xjt^ ■ (The expression for IC B I F m ^ used in our 
calculation is given in sect. 3.3.) 

Equations (32) and (34) imply the following form for 
the pion x amplitudes: 



X™(fc, fc ) = - E c n(W)V Nn (k) + V NN {k, ko) , 
n 

(35) 

X B ' T B {k,k u M',M) = -^c B (W,M)VM' n (k) 

K 

+V%P M {kM), (36) 



where V are the dressed vertices and V are the background 
parts of the amplitudes. 

The amplitudes involving the a-meson fulfill the same 
type of integral equations as the pion amplitudes with the 
kernels given in Appendix C. They assume the forms: 



X B Z(k,k M M,ri = -J2^ l (W,f J ,)V B I R (k)+V B f l ,(k,k li ) 
n 

(3 

X°j B T (k, h, M , M) = -*£c B (W, M)V^(fc) + Vf M (k, h) 
Xjt(Mm>m',m) = - 



n 



(37) 



(38) 



^n(W^)V^ K (k)+V^(k,k^. 

(39) 



n 



3.3 The Born approximation for the K matrix 

The Born approximation for the K matrix consists in ne- 
glecting the terms in (27), (28), (29), (32), and (34) in- 
volving the integrals. The K matrix is then constructed 
from the meson amplitudes (35)-(39) using (20) and (26), 
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and substituting the dressed vertices Vbti by the corre- 
sponding bare vertices Vbtz as well as V by JC. The ex- 
pressions for JC HH are derived in Appendix C; note that 
they involve only the on-shell amplitudes which are not in- 
fluenced by the approximation of separable kernels. They 
acquire the form 



\ - fB " 2M B ,. 
Jbb' 



B" 



2Ew[ 



+ M 2 „ -M 2 -ml 



and (in the Pll case) 



rB „ (h _ V- 2M B Vg, B (k li )V B , N (k 1 ) 

K-MpVtuW ~ 2^ 2Eu)ii + M 2 ; _ M 2 _ ^2 
= ^>m(^V' ^l) ' 

_ 2M B VS N (K)VS N {k ll ) 



(40) 



The symmetry in the background terms follows from the 
symmetry of the denominator in the u-channel, e.g. 
2E N uJ ti — Mfj — = 2E^lu — Mff — m 2 . Because we 
deal with s-wave scattering, the direct term (i.e. the first 
term in (37)) referring to the nucleon pole and the back- 
ground part for B = N almost cancel and can be dropped 
from the above sums. 

In the P33 case the decay into two correlated I = 
pions proceeds through the a A channel. As discussed in 
sect. 2.3 it is sensible to average over the A invariant mass 
such that the matrix elements depend only on the invari- 
ant mass of the two-pion system. The averaging is dis- 
cussed in the following. 



3.4 Averaging over invariant masses 

The averaging over the A invariant mass in the a A chan- 
nel implies that the matrix elements in the kernels assume 
the same form as (40) with TV replaced by A, and V^ B (k) 
by the averaged interaction matrix element defined as 



W-m„ 



UfiEV J M N +m 



(41) 

while the denominator in Kf" 7 contains lu^, k^ and 
evaluated by (15) in which Mm is replaced by the averaged 
invariant mass M. 

The averaged invariant masses M and p, are found 
by suitable smooth numerical approximations approach- 
ing the nominal hadron masses for large W, while remain- 
ing close to either Mm + or 2m n , for W slightly above 
the two pion threshold. 

The approximation of averaging over the invariant mass 
as in (41) can be applied to other matrix elements, as well 
as to the K matrix and the T matrix themselves. For the 
decay of a resonance B' into a pion and an unstable isobar 



B which in turn decays into the nucleon and the second 
pion, we introduce 



W(fci) 2 



kiW 

LUiE 



W-m„ 



dMwBiMfM^V^ik^ 



M N +m„ 



(42) 

where u>\ and E are those of (8) evaluated at M = M . 

Similarly, for the decay of a resonance B' through a 
baryon B and a cr-meson which in turn decays into two 
pions, the matrix element averaged over the meson invari- 
ant mass reads 

kW r w ~ MB 

VBB'Ck) 2 = ^ / d MWCT ( M ) 2 Ml V^ B ,(k,) 2 , 

where w = (W 2 - M% + p 2 )/2W, E = W - Q and k 2 = 
uj 2 +j2 2 . The averaging procedure turns the integral Hcitlcr 
equation (16) into a set of algebraic equations. Such an 
approximation does not influence the elastic channel since 
the matrix elements involving unstable hadrons in this 
channel always appear only under the integral; in inelastic 
channels this means that the M-dependent amplitudes are 
replaced by some averaged value. Identical averaging of 
amplitudes (42) is used in phenomenological analyses of 
ttN — > ttttN reactions proceeding through the unstable 
intermediate hadron (see e.g. [29]). 



3.5 Solving the integral equations in the 
approximation of separable kernels 

In this section we solve the set of coupled integral equa- 
tions (27), (28), (32), (34), beyond the Born approxima- 
tion for the K matrix. The solution for the vertices yields 
a considerable enhancement with respect to their bare val- 
ues while the solution for the coefficients c-jz involves a con- 
siderable mixing of different 'bare' resonances (denoted by 
1Z) appearing in our ansatz. The method yields simultane- 
ously the position of the resonance as well as the pertinent 
wave-function renormalization. Since the quark-u vertex 
is not as well determined as the quark- ir vertex, we treat 
the cr-meson vertices only in the Born approximation dis- 
cussed above. 

Inserting the ansaetze (35)-(39) into the set of coupled 
equations we obtain 



V Nn (k) = V NK (k) 
+ E 

B' 

^br(^) = V B n 



dk 



l K NN {k,k')V Nn {k') 
u' k + E N (k')-W 



dfc' M «' V - _ , (43) 



E B ,(k')-W 



BN 
M 



dk'^M 



(k,k')V Nn (k') 



E N (k') - W 



+ 



R' J 



, ^mm b , (^j k')V B ,\ (k') 
u>' k + E B .(k')-W 



dk 



■ (44) 



The background parts D (k, ko) and (k,ki) obey 
the integral equations of the type (43) with the non-homo- 
geneous terms JC NN (k, ko) and /C^ s (fc, fci), respectively, 
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while the T>f f N (k,k ) and T>f I P M (k 1 ki) satisfy the inte- 
gral equations of the type (44) with the non-homogeneous 
terms K.f t N (k,ko) and Kf { P M (k,k\), respectively. 

In appendix C we introduce several approximations 
which enable us to write the kernels in separable form: 



V B >N{k')V B >N{k) 



B' 

„ V B , N (k')V»? >B (k) 



(45) 



•BN tij 



— ^-m {k' i k) , 



(46) 



V#, B (fc')V#, B ,(fc) 



" K + e |„(M))K + e |;,(M'))' (47) 

where eg, = (M|, - M% - ml)/2E N and £§,(M) = 
(M|, - M 2 - ml) /IE, while /£ B are given by (65). Here 
M B < stands for the nominal (fixed) mass of the isobar B', 
while the invariant mass M pertinent to isobar B is a 
variable. Using the separable kernels we are able to solve 
the system exactly. This is important from the numerical 
point of view since it is now possible to control the princi- 
pal value integration over the poles of the kernel and thus 
avoid possible numerical instabilities. 

For the coefficients c|£ (here H denotes irN, nB, oB 
channels) we obtain a set of algebraic equations 



E A ™' (W)c%, (W, m H ) = b%(m H ) , 



(48) 



TV 



where 



ft' J 



'^MM B ,(k>kl)V B ,fc (k) 

Lu k + E B ,{k)-W 



= Kniki), 



(49) 



Here the sum over B' includes also the ground state, in the 
P33 case V a is replaced by Vak- The equalities between 
6 and V can be proved by iterating equations for V and 
V. 

Using the vector notation c H = [c|f, c^, , . . ] T and 
Vh = [Vhk, Vhk'i ■ • -] T : the solution of (48) can be writ- 
ten in the form c H = A~ 1 Vh ■ The zeros of A occur 
at the positions of the poles of the K matrix related to 
the resonance TZ; we denote the corresponding energies 
by M-r. The procedure to determine the coefficients cji is 



then the following: we first determine the zeros of the A- 
matrix determinant; by adjusting the energies of the bare 
states, Mjj, we can force the poles of the K matrix to ac- 
quire some desired values. (Note that in the case of several 
channels and strong background they do not coincide with 
the corresponding experimental values.) Diagonalizing the 
A matrix, UAU T = D, we write 



D = diag[A TC , \ n >, ■ ■ ■} 

= dmg[Z n (W)(W - Mi), Z n ,(W)(W - M R ,), 



(50) 



which defines the wave-function normalization Zr perti- 
nent to the resonance TZ. The solution can now be cast in 
the form 



c H = U t D _1 UVh • 



Finally, the resonant part of the x amplitudes appearing 
in the expression for the K matrix (e.g. (35)) takes the 
form 



X 



H'H 



-Vh,c m = -Vh.U 1 T> l XJVh 
_ S Vhk Z K (W)(W-M n ) 



Tl 



(51) 



n 



where 



Vhti — E u thvVhk' , c-n — 



TV 



Z n (W)(W -M n ) 



(52) 

The interpretation of (52) is that the resonant states TZ, 
1Z', ... are not eigenstates of H and therefore mix: 

<^ = E U TVJV$TV ■ 
TV 

The exception is the ground state - which by assumption 
is the eigenstate of H - for which the mixing of other res- 
onances because of (30) and (31) vanishes at the nucleon 
pole (W = Mn) and does not affect the ttNN coupling 
constant. 

Note that neglecting the off-diagonal terms A-rtv m 
(49) the expression for Zr in the vicinity of the resonance 
assumes the familiar form 



Zn{M n ) = 1 



W=M TL 



4 Results 

4.1 The parameters of the Cloudy Bag Model 

We analyze the capability of our approach in the frame- 
work of the Cloudy Bag Model (CBM) as one of the most 
popular examples of quark-pion dynamics. The Hamilto- 
nian of the model has the form (1) and (2) with 



1 



v{k) = 



k : 



,o 

■'mit 



3i(kR) 



2/ Vm^Ji wg, IT - 1 kR 
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assuming three quarks in the Is state. The parameter / 
corresponds to the pion decay constant /„., and w^it = 
2.043. It is a known drawback of the model that the ttNN 
coupling constant is underestimated, irrespectively of the 
bag radius, if / is fixed to the experimental value f v = 
93 MeV. We therefore adopt the conventional smaller value 
of / = 76 MeV which reproduces the irNN coupling con- 
stant. The free parameters of the model are the bag radius 
R and the energies of the bare quark states correspond- 
ing to the nucleon and the excited states. We have also 
considered alternative forms of the /c-dependence which 
avoid the typical oscillations due to the sharp cut-off at 
the bag surface but have found almost no change in the 
final results. We use the same bag radius for the excited 
states as for the ground state; as a consequence, the ma- 
trix elements of the quark-pion interaction between the 
quark configurations with different spatial structure are 
all proportional to v(k): 

{$B>\\V{k)\\cI> B ) 

3 

= r BB , r q v (k) ( J B , , T B , 1 1 ]T <r V 1 1 J B , T B ) , (53) 



Table 1. The model parameters used in sect. 4.2 (Born) and 
sect. 4.3 (Full) for the Pll and P33 partial waves. Mb are the 
positions of the K-matrix poles, Mb* corresponds to iV(1710). 
The parameter r B B' is defined in (53) and f^ BB i refers to the 
values used in the kernels in the full calculation. 





Born 


Full 


Parameter 


Pll P33 Pll 


P33 


R 


0.83 fm 




M A 


1232 MeV 




Mr 


1520 MeV 




Ma- 


1780 MeV 




Mr, 


1870 MeV 






450 MeV 




r a 


550 MeV 






0.96 


0.99 


TNR 


1.68 0.75 


1.00 


tna 


0.83 1.39 


1.00 


TAR 


0.80 


1.30 


fnNA 


1.30 


1.12 


TnAA 


1.00 


1.25 


r-KNN 


1.00 


1.00 



1 if both B and B' are in the (Is) 3 configu- 
= r w for the transition between the (ls) 2 (2s) 1 

- r 2 , if 



where r q 
ration, r q 

configuration and the ground state, and r q = 
both B and B' are in the (Is) 2 (2s) 1 configuration. Here 

1 I _ i \ 1 1/2 

-"mitI^mit l ) 



1 

71 



U) 



\:It( w MIT 1) 



with wJj IT = 5.396. The parameter r BB > in (53) allows us 
to tune the chosen coupling constant relative to its ST/ (6) 
value. The value r^N — 1 is fixed by our choice of /. We 
assume r BB i = r B i B . 

In the Pll case the sum over TZ in (22), (23) and (24) 
includes beside the nucleon, the Roper iV(1440) and the 
V(1710), and in the P33 case, the 4(1232), 4(1600) and 
4(1920). We do not include further intermediate states 
and channels since our present goal is to find a pattern 
common to the low-lying Roper-like resonances. This lim- 
its the validity of our approach to energies below — 1700 MeVthe energies close to the 4(1600) resonance. This is most 



4.2 The Born approximation 

We first analyze the P33 partial amplitudes. The results 
turn out to be almost insensitive to the value of the bag 
radius so the only parameters to adjust are the positions 
of the resonances Ma and Ma* and the relative coupling 
strengths r^A and r^A* = t^atnr- (Note that our value 
for Ma* is the position of the K-matrix pole and should 
not be identified with the nominal value of the resonance 
invariant mass.) 

In order to investigate the importance of different de- 
grees of freedom we include in the first step only the irN 
and the 7r4 channels without the background. In this case 
the model reproduces the amplitudes at lower energies 
provided we take a value for r^A which is substantially 
larger than unity (Fig. 1). When the background is in- 
cluded, the agreement considerably improves except for 



The approach can be extended in a straightforward way 
by including higher intermediate states as well as other 
channels. 

In sect. 4.2 we first discuss the results of the Born 
approximation for the K matrix, and in sect. 4.3 the results 
when the integral equation for the K matrix is solved. The 
parameters in the two cases are displayed in Table 1. The 
parameters Ma, fn a and F a are kept fixed: Ma at the 
experimental position of the pole of the K matrix, while 
from the recent analysis of Leutwyler [57] we take m a = 
450 MeV and R a = 550 MeV. The values for M R , M A * 
and Mji- , and for G a are free in the Born approximation 
as well as in the full calculation. The parameters r BB > 
defined in (53) are free in the Born approximation; in the 
full calculation they are kept at the values predicted by 
the quark model except for the value of the nAR coupling. 
The 7Vb b correspond to the averaged values of the dressed 
vertices and are explained in sect. 4.3. 



notably seen in Fig. 2 in which the inelasticity exhibits a 
qualitatively different behavior compared to the case with 
no background, and becomes consistent with the experi- 
mental data up to W ~ 1700 MeV. A small kink around 
1700 MeV is an indication of the 4(1600) resonance. 

Including the a A and the ttR channels gives an al- 
most perfect fit to the experimental amplitudes also in 
the vicinity of the 4(1600) resonance, washing out almost 
completely the signature of that resonance in the phase 
shift (Fig. 1). 

Turning to the case of the Pll partial wave we note 
that including only the 7rV and it A channels without the 
background the Born approximation fails to reproduce the 
amplitudes determined in the partial wave analysis even 
if we considerably increase the values of some coupling 
constants and take = 2.20 and tar, = 1-55 (Figs. 1 
and 4) . Adding the background yields the correct behavior 
of the amplitudes below the two-pion threshold. In order 
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5 (P33) 



nel reproduces the threshold behavior and considerably 
improves the agreement at higher energies. 



tt/2 




1200 



1400 



1600 



1800 
W [ MeV ] 



Fig. 1. The P33 (top panel) and Pll (bottom panel) phase 
shifts in various approximations. Born approximation: nN and 
it A channels, resonant terms only (dotted lines), with back- 
ground (dashed lines), adding the a channel (thin solid lines). 
Full calculation: thick solid lines. The data points in this and 
subsequent figures arc from the SAID ttN — ► ttN partial-wave 
analysis [32] unless noted otherwise. The model parameters are 
given in Table 1. 





1200 



1400 1600 1800 

W [ MeV ] 



Fig. 2. Inelasticity in the P33 (top panel) and the Pll wave 
(bottom panel). Notation is as in Fig. 1. 



to reproduce the amplitudes in the region of the Roper 
resonance we have to keep the large value for tar, still, 
the approximation does not reproduce the rapid rise of 
the inelasticity just above the two-pion threshold nor the 
property that it remains close to unity even well above the 
Roper resonant energy (Fig. 2). Including the aN chan- 




1600 1800 
W [ MeV 1 

Fig. 3. The real and imaginary parts of the T matrix for the 
P33 partial amplitudes. Notation as in Fig. 1. 



0.6 
0.4 
0.2 
-0 
-0.2 
-0.4 



1 

■ p r 


1 1 

ImT : 








\\ T^>x 







1200 



1400 



1600 1800 
W [ MeV ] 



Fig. 4. The real and imaginary parts of the T matrix for the 
Pll partial amplitudes. Notation as in Fig. 1. 



Let us stress that in order to reproduce the widths of 
the AT(1440) and Z\(1232) resonances, the respective irNR 
and 7riV A coupling constants have to assume considerably 
larger values than the corresponding bare-quark values. 
The choice of the model parameters tbb 1 (53) that drive 
these constants turns out to be quite different for the P33 
and the Pll case. A possible solution to these inconsis- 
tencies is presented in the next section. 



4.3 Results of the full calculation 

The results presented here have been obtained using the 
same set of parameters as in the Born approximation for 
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the bag radius, the mass and the width of the cr-meson, 
and the position of the K-matrix poles (see Table 1). Sim- 
ilar results are obtained for 0.75 fm < R < 1.0 fm as well 
as for 400 MeV < m CT < 550 MeV provided the coupling 
constants are slightly readjusted. We keep the same set 
of parameters for both partial waves; we allow only for a 
slight deviation in the parameters entering the kernels. 

We first investigate the relation between the bare ma- 
trix elements of the baryon- meson interaction (53) and the 
corresponding dressed values Vbb> which are solutions of 
the integral equations (43) and (44). In the case of the 
elastic channel we introduce the ratio 



T-nNB 



VjVB(fcp) 

V/vs(fco) 



(54) 



measuring the renormalization of the bare vertex. In the 
case of the iV(1440) and the Z\(1232) we find that this ra- 
tio exhibits a relatively strong energy dependence (Fig. 5) 
and yields a substantial enhancement of the bare coupling 
constant over a broad energy range. The enhancement is 
consistent with the value of the corresponding coupling 
constant used in our analysis in the Born approximation. 

The formulas (43) and (44) represent a system of cou- 
pled non-linear integral equations for the dressed vertices 
since the kernels (45) - (47) themselves contain these ver- 
tices. We have not attempted to solve the system exactly 
but have substituted the dressed vertices appearing in the 
kernels by VsB'ik) — r^BB'VBB'ik) with suitably chosen 
values for t^bb 1 ■ We have adjusted these values by averag- 
ing the corresponding solutions for the dressed vertices in 
the relevant energy range, allowing for small variations to 
obtain better overall fits. This approximation is justified 
because the contribution from the integrals in (43) and 
(44) turns out to be less important compared to the lead- 
ing term. Also, many of the dressed vertices only negligibly 
influence the result such that we can put the correspond- 
ing fnBB' = 1. In Table 1 only a few relevant cases are 
listed; in all other cases t^bb 1 = 1- 

In the Pll wave the calculated amplitudes closely fol- 
low the experimental values in the energy range from the 
threshold up to W ~ 1700 MeV (Figs. 4 and 2). In particu- 
lar, we should stress the excellent agreement of the ampli- 
tudes slightly above the 2tt threshold where the inelasticity 
is dominated by the aN channel. Compared to the results 
of the Born approximation we notice an improvement at 
very low energies as well as at those above the resonance, 
which is a consequence of the energy-dependent dressing. 

Above W ~ 1700 MeV the imaginary part of the am- 
plitude as well as the inelasticity drop rapidly. The simple 
model involving only the Roper-like bound state and the 
nucleon breaks down here and the effects of other reso- 
nances as well as other channels such as the pN and r]N 
may become more important. To investigate this point we 
have included in our calculation the next Pll excitation, 
the JV(1710), and assumed that it couples only to the aN 
channel, representative of a generic irir decay with ~ 80 % 
of the value of G a . (The experimental branching fraction 
for jV(1710) -> Nittt is (40-90) %.) A better agreement at 
higher energies is obtained, although a similar effect can 
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Fig. 5. Ratios (54) of the 7riVii (solid line) and the ttNA 
(dashed line) transition amplitudes to the respective bare 
quark values as a function of the invariant mass W. Straight 
lines: corresponding constant values used in the Born approx- 
imation. Model parameters as in Table 1. 
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Fig. 6. Inelastic contributions to the imaginary part of the T 
matrix for the Pll partial wave. Full lines and filled circles: con- 
tribution of AT(1440) — > NIivkYsZwuvc- Dashed lines and empty 
circles: contribution of iV(1440) — ► An. Thin lines represent 
the Born approximation. The data are from [29]. Model pa- 
rameters as in Table 1. 



be achieved by decreasing the bag radius (or increasing 
the cut-off parameter in general) which makes a model- 
independent analysis less reliable. 

In Fig. 6 we display the imaginary parts of the off- 
diagonal T- matrix elements corresponding to the processes 
irN -> irA and irN -> aN calculated with the JV(1710) in- 
cluded. The amplitudes are averaged over the range of un- 
stable hadron invariant masses as described in sect. 3.4. In 
comparison to the values extracted from a recent partial- 
wave analysis of irN — > ttN and irN — > ttttN [29], our 
model overestimates the aN decay probability and un- 
derestimates the ttA channel. The disagreement is a con- 
sequence of the destructive interference between the nu- 
cleon and the bare Roper in this channel at higher W. 
In the Born approximation, a better agreement in both 
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inelastic channels is obtained. Had we departed from our 
standard set of parameters chosen to yield a good overall 
description of various amplitudes in both partial waves, 
or by including further channels in the integral equations, 
an improvement can be achieved in the full calculation as 
well. (Note also the large systematic scatter of data.) 

The full calculation for the P33 partial wave does not 
significantly improve the results of the Born approxima- 
tion which are anyway satisfactory throughout the energy 
range from the threshold to W <~ 1700 MeV. It does, how- 
ever, explain the strong enhancement of the bare quark- 
model ttNA coupling of the Born approximation. Above 
W <~ 1800 MeV the agreement is lost even if we include 
another bare resonant state corresponding to the Z\(1910), 
indicating the need to include further channels not present 
in our analysis. 
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Fig. 7. The real and imaginary parts of the T matrix and 
the inelasticity for the Pll partial wave for different bag radii: 
R = 1 fm (doted lines), R = 0.9 fm (dashed lines), R = 0.83 fm 
(solid lines), and R = 0.77 fm (dashed-dotted lines). 

In order to test the sensitivity of the results we have 
varied the parameters appearing in Table 1. Varying the 
bag radius and the parameters r^B we find that in the 
P33 case the results are almost independent of the bag ra- 
dius while a few percent change in r^B 's may improve the 
agreement. In the Pll case, the results are more sensitive 
to the choice of the bag radius (Fig. 7). For R ~ 0.8 fm the 
agreement at larger energies improves; for R < 0.75 fm the 
real part of the amplitude comes too close to zero and the 
phase shift loses its characteristic resonant behavior. The 
optimal value of r^R increases by ~ 5 % for R — 1.0 fm 
and decreases by <~ 10 % for R = 0.77 fm. 

In Table 1 we list the positions of poles of the K matrix 
which are free parameters in our approach. Alternatively, 
we could have used the masses of the bare quark states en- 
tering our ansaetze for the PV states (22), (23) and (26); 
these masses are related to the pole values through the 
procedure described in sect. 3.5. From the values of the 
bare masses we may obtain an indication of how sensible 
is the comparison of baryon masses calculated in various 



quark models (with no mesons included) and the experi- 
mental positions of the corresponding resonances. We are 
not able to give a conclusive answer since the result is 
strongly sensitive on some quantities (in particular the 
bag radius and the irRR coupling constant) which only 
weakly influence the behaviour of the scattering ampli- 
tudes and are therefore not well fixed in our calculation. 
Nevertheless, we can conclude that the self-energy correc- 
tions and the effects of resonance mixing most strongly 
affect the ground state, lowering its mass compared to 
the bare mass by 310 MeV at i?b ag = 0.83 fm (or by 
180 MeV at i?bag = 1-0 fm). The bare delta-nucleon mass 
difference turns out to be smaller than the experimen- 
tal one by about 90 MeV; the remaining difference can 
be attributed to the chromomagnctic interaction. In the 
case of the Roper resonance the bare mass splitting is 
smaller by 120 MeV at R hag = 0.83 fm (or by 60 MeV 
at i?ba g = 1.0 fm). Since our value for the K matrix pole 
is typically 100 MeV higher than the experimental value, 
the resulting bare Roper-nucleon mass splitting remains 
in the ballpark of admissible values. A similar conclusion 
holds for the Z\(1600) where the bare mass splitting is typ- 
ically 200 MeV lower than the value deduced from the K 
matrix pole. 

5 Conclusions and perspectives 

We have developed a general method to incorporate the 
excited baryons represented as quasi-bound quark-model 
states into a dynamical framework with correct unitarity 
and symmetry requirements and proper boundary con- 
ditions. To illustrate the key points, we have used the 
Cloudy Bag Model, although the method is applicable to 
more sophisticated models. In fact, the only information 
needed as input to our computational scheme are the ma- 
trix elements of the meson interaction between the quark 
states, i.e. (4>(k)\V\<l>(k = 0)). In a more ambitious com- 
putational scheme, the method can be extended to allow 
for the readjustment of the intrinsic three-quark wave- 
function to the scattering boundary conditions by using 
the Kohn variational principle. 

We have shown that an intricate interplay of the wN, 
irA and <jN degrees of freedom governs the elastic and 
inelastic pion-nucleon scattering in the energy range from 
the threshold up to W ~ 1700 MeV. The model explains 
at least qualitatively the behavior of the amplitudes in the 
vicinity of the Z\(1600) which can be regarded as the P33- 
wavc counterpart of the iV(1440) in the Pll wave. We have 
studied in a systematic way the role of the background 
processes which turn out to be important to qualitatively 
reproduce the experimental data throughout the first and 
the second resonance region. 

We have described the correlated two-pion decay in 
the relative s-wave by the cr-meson. In spite of this purely 
phcnomenological approach, our results show that this de- 
gree of freedom is crucial to explain certain features of the 
scattering amplitudes, in particular the inelasticity for the 
Pll partial wave just above the two-pion threshold which 
rapidly rises from zero to unity and remains close to this 
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value in a broad energy range. Because of the s-wave na- 
ture of the cr-meson coupling compared to the p-wave cou- 
pling of the pion to the A, the two-pion decay is dominated 
by the aN channel in the energy region slightly above the 
two-pion threshold. This offers the possibility to determine 
the strength of the coupling in a model-independent way, 
which is generally not the case at higher energies where 
the same level of agreement can be obtained in a broad 
range of parameters. Our calculation indicates that a low 
mass of the cr-meson (m a m 450 MeV) is preferable over 
the larger mass typically obtained in partial-wave analy- 
ses. 

We have designed a framework to numerically solve 
the integral equation for the K matrix by approximat- 
ing the kernels with a separable form which preserves the 
symmetries of the matrix and thus ensures unitarity. The 
important outcome of this calculation is a substantial in- 
crease of the quark-model pion-baryon couplings explain- 
ing the large width of the Roper and possibly other reso- 
nances. Such a strong enhancement could — at least par- 
tially — accommodate relatively weak pion-baryon cou- 
pling strengths predicted in constituent quark models. This 
shows, on the one hand, that one can explain the nature of 
the Roper resonances without invoking "exotic" degrees of 
freedom mentioned in the Introduction and, on the other 
hand, establishes one of the benchmarks for an assess- 
ment of the underlying quark models. Still, the scattering 
analysis alone should not be expected to provide a defini- 
tive selection criterion, because a consistent description, 
as shown in the example above, can be achieved within 
a relatively broad range of parameters. A future analysis 
devoted to pion electro-production could provide a more 
complete set of criteria. The application of the method to 
calculate the electro-production amplitudes will be treated 
in a separate paper. 



A Evaluation of one- and two-pion matrix 
elements 

We derive some expressions for the matrix elements of the 
pion field between the cigcnstatcs of the Hamiltonian with 
the pion part given by (1). If & A is an eigenstate then 

(w fc + H - E A )a mt {k)\* A ) = -Vl t (k)\* A ) , (55) 



(ui k + u' k + H - E A )a mt {k)a m > t >{k')\\P A } 

= (yUk)a m >Ak') + Vl, t ,(k')a mt (k)) Wa) , (56) 

Multiplying (55) by an eigenstate of (1) we obtain 

($ B (k)\a mt (k)\V A ) = v5{uj k + E B (k) - E A ) 

(*B{k)\VUWA) 



{uj k + E B {k)-E A ) 



(57) 



where v is an arbitrary constant; it is used to determine 
the normalization of the channel states (19), (23) and (24). 



Taking for \P A either (19) or (23), multiplying (56) by 
an eigenstate with momentum k + k' , and neglecting the 
terms with two or more pions, we obtain 

(w fc + u' k + E B (k + k') - W)$ B {k + k')\a mt {k)\y B ,{k')) 
= -$ B {k + k')\vUk)WB>{k')). (58) 

Note that when B' is on-shell, E B ' = W — uj' k , and (58) 
reduces to (57). 



B Derivation of the mass distribution function 

We derive here the expression for the mass distribution 
function w B (M) in (21) for the case of the ttA channel 
which dominates the two-pion decay through the inter- 
mediate baryon. In this case we can assume that the in- 
termediate A decays only into a pion and the nucleon. 
(This assumption is justified since experimentally the elas- 
tic channel remains the dominant process also well above 
the two-pion threshold.) Then the K matrix becomes a 
scalar (denoted as Ka) and the orthonormalized state in 
(10) assumes the form 



Wa{M)) 



1 



1 



w a {M)\$a{M)) . 



(59) 



where M is the invariant mass of the intermediate A, while 
u> 2 and fc 2 are the energy and momentum of the second 
pion. From (35) we see that close to the resonance the pion 
amplitude behaves as x A N (k, k 2 ) 



-c^(M)V N A(k) and 



from (36) we find Xa (k,k 2 ,M B ) 



(M) V 



EA CO' 



(M - MaY 1 - The state & A {M) introduced in 



with 

(59) can then be written in the form 



Wa{M))=za{ \$a) 



- dk 



V NA (k) 



uj k + E N (k) 
-£/d/ 



M 



v£f(fc) 



u k + E B (k) - M 



[at(fc)|<W 



[at(fc)|# B )]i 



(60) 



Note that in the weak-coupling limit the form (60) cor- 
responds to the usual perturbative expression for the A 
state, with tf^v and &a replaced by the corresponding bare 
quark states. It is dominated by the bare-quark state <Pa ■ 
This is another reason for choosing the particular form of 
factorization in (59). 

From (51) and (26) we have (K A = Kp N ) 



K A (M) = TT 



u) 2 M N V N A{k2)' 2 



k 2 M Z A (M)(M A - M) 



(61) 
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for M w M A , where the terms denoted by . . . vanish at 
the resonance; Z A {W) is defined in (50). From (59) it then 
follows 



we substitute the momenta of the intermediate baryon by 
a suitably chosen average over k + fc' denoted as k. The 
kernel then takes the form: 



w A (M) 



1 



fc 9 M 



1 



Z A y/1 + K\ 7T V UlM N V NA {k 2 ) ' 



At this point we can use the expression for the K ma- 
trix obtained numerically or use a suitable parameteriza- 
tion either of the computed form or of the experimental 
data. The simplest choice is 



C 



M A -M' 



c 



ujjM N V NA (kj) 2 
k£M A Z A (M A ) 



(62) 



where u£ = (M\ - M% + ml)/2M A and the residue C 
(corresponding to \r) is assumed to be VT-independent; 
the second equality follows from (61). From (21) and (60) 
we have 



w A {M) = 




z Ax /Z A (M A ) 



For sufficiently small values of C, w A (M) in (63) is strongly 
peaked around M A and has a unit integral provided that 
hence 



z A \J Z A (M A ) 



1 



w A {Mf -> S(M - M A ) , C -» . 



(64) 



In this limit we obtain the usual wave-function normaliza- 
tion of (60), i.e. z A = Z A (M A )-V 2 . 

In a more precise calculation we use a better approxi- 
mation for the K matrix, K a = C/(M A — M) + D, which 
turns out to give a very good approximation for M from 
the threshold to values well above the resonance. Note, 
however, that the parameter D should not be identified 
with the background contribution since in the approxi- 
mate formula (62) the coefficient C is kept fixed while 
in (61) all terms exhibit strong k 2 (or, cquivalcntly, M) 
dependence. The coefficient D is chosen such as to com- 
pensate this dependence. In our calculation we use C = 
55 MeV and for D a constant value —0.41 below M <~ 
1400 MeV and, above it, a value that smoothly approaches 
zero. 



C Evaluating the kernels 

To evaluate the kernel (33) in the integral equation for 
Xjt amplitudes we insert a complete set of states 



1 



\*n){*n\+z2 \*b(M))(* b (M)\ 



where the sum implies also the integral over invariant 
masses and momenta. Next we use the (adjoint of) (57) 
and (58). Following Ericson and Weise (see [49], sect. 2.5.3.) 



(, US _ f N (*B\\V(k)\\* N }(* N \\V(k'W N ) 
[K, K ) — J NB ■ 



LU k +Lu' k +E N {k)-W 



■£/. 



B 

NB 



dM 



, (^nnfc)ii^(MQ) 

cj k +cj' k + E'(k)-W 
x(V N \\V(k')0 B ,(M')), 



where 



Jab 



V(2J A + 1)(2J B + 1)(2T A + 1)(2T B + 1) 

x W(1J a JbU J c , J)W(1T a T b 1; T c , T) . (65) 



The matrix elements are those entering the expression 
for the x's in (35)-(36), {&N\\V(k')\\& B ,(M')) = w B ,{M') 

V B , N (k') and (§ B (M)\\V\\$ B ,(M')) - w B ,(M')V^ B . We 
can now approximately perform the integration over M' 
assuming (64): 



V^ B (k)V NN (k>) 
lNB u k +u/ k + E N (k)-W 



+ £ Inb 



V^ B (k)V B , N (k') 
uj k +uj' k + E B ,(k)-W 



+ ... 



In order to build up a feasible computational scheme we 
approximate the kernel with a separable expression, a rel- 
ativistic extension of the approximation used in Ref . [23] : 

1 



w fc + w' k + E B > (fc) - W 

lj + wi + E B > (fc) - W 



E B ,(k) - E N (k ))(co k + E B ,{k) - Eih)) 



Here k\ , uj\ are the on-shell pion momentum and energy in 
the ttB channel satisfying W = u)\+E(ki) = wo + £^v(fco)- 
The approximation on the RHS coincides with the exact 
expression on the LHS when either of the two pions is on- 
shell, i.e when either ui k = ujq or w' k = u>\. (This can be 
easily seen by writing W — ojq + £jv(fco) on the LHS and 
W = lo\ + E(k\) on the RHS in the first case and vice 
versa in the second one.) When both pions are on-shell, 
the denominator can be cast in the form 

1 E B >(k)+E N (ko)-e>i 



u 1 + E B '(k)- E N (k ) 



Ml, 



2u! 1 E N (k ) - - ml 
(66) 

We have assumed (k ■ ki) — which is essentially the 
same approximation suggested in [49] and the denomi- 
nator acquires the characteristic u-channel form. Since we 
describe the resonance in terms of the wave function rather 
than bispinors, our numerator differs from the correct rel- 
ativistic expression in the u-channel (see e.g. [49]). Since 
our expression is anyway approximate, we replace the fac- 
tor E B '(k) + E N (k ) -u>i = E B '(k) + E B {ki) - oj in our 
numerator by the correct relativistic expression 2M B >. 
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In the general case we find: 



V B"B 



uj k +u' k + E B »(k) - W 



where the sum includes also the ground states. We now 
approximate 



uj k +u' k + Eb" (k) 



W 



E B „{k)-W 



(w' k +E B „{k)-E{k 1 )) (w fc + £ B " (k) - E'(k[ )) 

Here hi, UJ\, k[, u>[ are the on-shell pion momenta and 
energies satisfying W = <jj\+E(k\) = lo[+E' (k[), E(ki) = 



y^M' 2 + k[ 2 . The approximation on 



VM 2 + k 2 , E'(k[ 

the RHS of (67) coincides with the exact expression on the 
LHS of (67) when either pion is on the mass shell, i.e. when 
cither io k = u>[ or uj' k = uj\. When on-shell, the kernels are 
proportional to the u-channel background elements of the 
K matrix. We now use the same approximation as in (66) 
for k as well as for the numerator; the final forms are given 
by (45)-(47). 

In the case of the channels involving the cr-meson wc 
find 



ft.k 



■to k + E B >(k)-W 



W 7 

B' M fe 



+ u^ k + E B i{k)-W 



where in the Pll case the sum over B 1 includes only 
J = T = i isobars; in the Pll case the N in the above 
expression is replaced by A while the sum over B' includes 
only J = T = | isobars. 



References 

1. T. Melde, W. Plessas, R. F. Wagenbrunn, Phys. Rev. C 72 
(2005) 015207; T. Melde, W. Plessas, R. F. Wagenbrunn, 
Phys. Rev. C 74 (2006) 069901(E). 

2. Y. B. Dong, A. Faessler, K. Shimizu, Nucl. Phys. A 705 
(2002) 153. 

3. R. Koniuk, N. Isgur, Phys. Rev. D 21 (1980) 1868; R. Ko- 
niuk, Nucl. Phys. B 195 (1982) 452. 

4. C. P. Forsyth, R. E. Cutkosky, Z. Phys. C 18 (1983) 219. 

5. N. Isgur, G. Karl, Phys. Rev. D 18 (1978) 4187; Phys. Rev. 
D 19 (1979) 2653. 

6. S. Capstick, W. Roberts, Phys. Rev. D 47 (1993) 1994. 

7. L. Theussl, R. F. Wagenbrunn, B. Desplanques, 
W. Plessas, Eur. Phys. J. A 12 (2001) 91. 

8. F. Stancu, P. Stassart, Phys. Rev. D 38 (1988) 233. 

9. T. Burch el al, Nucl. Phys. A 755 (2005) 481. 

10. D. Guadagnoli, M. Papinutto, S. Simula, Nucl. Phys. A 
755 (2005) 485; Phys. Lett. B 604 (2004) 74. 

11. N. Mathur et al, Phys. Lett. B 605 (2005) 137. 



12. S. Sasaki, T. Blum, S. Ohta, Phys. Rev. D 65 (2002) 
074503; K. Sasaki, S. Sasaki, T. Hatsuda, Phys. Lett. B 
623 (2005) 208. 

13. M. Dillig, M. Schott, Phys. Rev. C 75 (2007) 067001. 

14. K. Khosonthongkee, N. Supanam, Y. Yan, Th. Gutsche, 

A. Faessler, Nucl. Phys. A 790 (2007) 518c. 

15. B. Julia-Di'az, D. O. Riska, Nucl. Phys. A 780 (2006) 175. 

16. P. Alberto, M. Fiolhais, B. Golli, and J. Marques, Phys. 
Lett. B 523, 273 (2001) 

17. J. A. Elsey, I. R. Afnan, Phys. Rev. D 40 (1989) 2353; 

B. C. Pearce, I. R. Afnan, Phys. Rev. C 40 (1989) 220. 

18. E. D. Cooper, B. K. Jennings, Phys. Rev. D 33 (1986) 
1509. 

19. E. A. Veit, B. K. Jennings, A. W. Thomas, Phys. Rev. D 
33 (1986) 1859. 

20. G. Kalbermann, J. M. Eisenberg, Phys. Rev. D 29 (1984) 
517. 

21. A. Suzuki, Y. Nogami, N. Ohtsuka, Nucl. Phys. A 395 
(1983) 301. 

22. A. S. Rinat, Nucl. Phys. A 372 (1982) 341. 

23. A. W. Thomas, S. Theberge, G. A. Miller, Phys. Rev. D 
24 (1981) 216; S. Theberge, A. W. Thomas, G. A. Miller, 
Phys. Rev. D 22 (1980) 2838; erratum Phys. Rev. D 23 
(1981) 2106. 

24. Z. Li, V. Burkert, Z. Li, Phys. Rev. D 46 (1992) 70; 

C. E. Carlson, N. C. Mukhopadhyay, Phys. Rev. Lett. 67 
(1991) 3745. 

25. W. Broniowski, T. D. Cohen, M. K. Banerjee, Phys. Lett. 
B 187 (1987) 229. 

26. E. Hernandez, E. Oset, M. J. Vicente Vacas, Phys. Rev. C 
66 (2002) 065201. 

27. N. Fettes, U.-G. Meifiner, Nucl. Phys. A 693 (2001) 693. 

28. T. S. Jensen, A. F. Miranda, Phys. Rev. C 55 (1997) 1039. 

29. D. M. Manley, E. M. Saleski, Phys Rev. D 45 (1992) 4002 
and references therein, and private communication. 

30. R. E. Cutkosky, S. Wang, Phys. Rev. D 42 (1990) 235. 

31. R. S. Longacre, J. Dolbeau, Nucl. Phys. B 122 (1977) 493. 

32. R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, R. L. Work- 
man, Phys. Rev. C 74 (2006) 045205. 

33. T. P. Vrana, S. A. Dytman, T.-S. H. Lee, Phys. Rep. 328 
(2000) 181. 

34. S. Ceci, A. Svarc, B. Zauner, Phys. Rev. Lett. 97 (2006) 
062002. 

35. A. Matsuyama, T. Sato, T.-S. H. Lee, Phys. Rep. 439 
(2007) 193; 

36. B. Julia-Diaz, T.-S. H. Lee, A. Matsuyama, T. Sato, Phys. 
Rev. C 76 (2007) 065201. 

37. H. Kamano, M. Arima, Phys. Rev. C 73 (2006) 055203; 
H. Kamano, M. Morishita, M. Arima, Phys. Rev. C 71 
(2005) 045201; H. Kamano, M. Arima, Phys. Rev. C 69 
(2004) 025206. 

38. G. Penner, U. Mosel, Phys. Rev. C 66 (2002) 055211. 

39. O. Krehl, C. Hanhart, S. Krewald, J. Speth, Phys. Rev. C 
62 (2000) 025207. 

40. H. P. Morsch, P. Zupranski, Phys. Rev. C 61 (1999) 
024002. 

41. A. B. Gridnev, N. G. Kozlenko, Eur. Phys. J. A 4 (1999) 
187. 

42. T. Feuster, U. Mosel, Phys. Rev. C 59 (1999) 460. 

43. C. Schiitz, J. Haidenbauer, J. Speth, J. W. Durso, Phys. 
Rev. C 57 (1998) 1464. 

44. F. Gross, Y. Surya, Phys. Rev. C 47 (1993) 703. 



B. Golli, S. Sirca: Roper resonances in chiral quark models 



17 



45. P. Alberto, L. Amoreira, M. Fiolhais, B. Golli, S. Sirca, 
Eur. Phys. J. A 26 (2005) 99. 

46. S. Schneider, S. Krewald, Ulf-G. Meifiner, Eur. Phys. J. A 
28 (2006) 107. 

47. R. G. Newton, Scattering Theory of Waves and Particles, 
Dover Publications, New York 1982. 

48. G. F. Chew and F. E. Low, Phys. Rev. 101, 1570 (1956). 

49. T. Ericson and W. Weise, Pions and Nuclei, Clarendon 
Press, Oxford 1988. 

50. M. C. Birse and M. K. Banerjee, Phys. Lett. B 136 (1984) 
284; Phys. Rev. D 31 (1985) 118. 

51. S. Kahana, G. Ripka, and V. Soni, Nucl. Phys. A 415 

(1984) 351. 

52. W. Broniowski and M. K. Banerjee, Phys. Lett. B 158 

(1985) 335; Phys. Rev. D 34 (1986) 849. 

53. C. V. Christov, A. Blotz, H.-C. Kim, P. V. Pobylitsa, T. 
Watabe, Th. Meissner, E. Ruiz Arriola, and K. Goekc, 
Prog. Part. Nucl. Phys. 37 (1996) 1. 

54. G. Ripka, Quarks Bound by Chiral Fields, Clarendon 
Press, Oxford 1997. 

55. W. Broniowski and B. Golli, Nucl. Phys. A 741 (2003) 
575. 

56. S. Eidelman et al. (Particle Data Group), Phys. Lett. B 
592, 1 (2004). 

57. H. Leutwyler, Int. J. Mod. Phys. A 22 (2007) 257. 



